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ABSTRACT 

We investigate the relative distribution of the gaseous contents of the Universe 
(as traced by a sample of Lyman a (Ly-a ) absorbers), and the luminous baryonic 
matter (as traced by a redshift survey of galaxies in the same volume searched for 
Ly-a absorbers), along 16 lines-of-sight (LOS) between redshifts and 1. Our galaxy 
redshift survey was made with the Multi-Object Spectrograph (MOS) on Canada- 
France-Hawaii Telescope (CFHT) and, when combined with galaxies from the litera- 
ture in the same LOS, gives us a galaxy sample of 636 objects. By combining this with 
an absorption line sample of 406 absorbing systems drawn from published works, we 
are able to study the relationship between gas and galaxies over the latter half of the 
age of the Universe. A correlation between absorbers and galaxies is detected out to 
separation of 1.5 Mpc. This correlation is weaker than the galaxy-galaxy correlation. 
There is also some evidence that the absorbing systems seen in CIV are more closely 
related to galaxies, although this correlation could be with column density rather 
than metallicity. The above results are all consistent with the absorbing gas and the 
galaxies co-existing in dark matter filaments and knots as predicted by current mod- 
els where the column density of the absorbing gas is correlated with the underlying 
matter density. 

Key words: galaxies - intergalactic medium, galaxies - quasars: absorption lines, 
galaxies - galaxies: haloes. 



1 INTRODUCTION 

This paper is part of our efforts to provide constraints and 
measurements of the relationship between the distribution 
of the gaseous contents of the Universe as traced by neutral 
Hydrogen and the luminous baryonic mater as traced by 
galaxies. 

The manner in which gas collapses gravitationally into 
dark matter potential wells to form stars (and hence galax- 
ies), and the way in which these stars then affect the gas is 
a topic of great interest at present. This is being pursued 
both observationally and theoretically at redshifts from 7 
to zero (i.e. when the universe was approximately 5% of its 
current age, to the present day). In this paper we investi- 
gate this process observationally during the second half of 
the evolution of the universe. 
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There is a continuing debate about the relationship be- 
tween low redshift Lyman a (Ly-a) absorbers and galaxies, 
where here 'low redshift' is taken to mean redshifts less than 
1. A simplified (strawman) version of the two sides is (a) 
that all low redshift Ly-a absorbers are part of physically- 
distinct luminous-galaxy halos, or (b) that they are all part 
of the filamentary structure seen in recent SPH/Mesh struc- 
ture formation models, and are only related to galaxies by 
the fact that both are following the underlying dark mat- 
ter distribution. In practice, almost all authors acknowledge 
that the universe includes a mix of the above two popu- 
lations (and indeed others), and the debate is more about 
which population dominates a particular set of observations. 

Both of these positions have been vigorously defended in 
the literature. The reference list below includes all relevant 
papers listed on the ADS abstract server from the time of 
the first available high quality UV spectroscopy from HST 
to 15 August 2005. 

Paper s supporting (a) include ( in chronological order) : 
|MoI J1994): lLanzetta et all il995T) : lLanzetta et al] |l99rji : 
Barcons et alJ dl99Slh iLinderl il998l> : IChen et*afl !l99cl) : 
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Ortiz-Gil et alJ Il999h: Hinder! feOOoh: IChen et all l200lft: 



Chen et all fcOOll): ISteidel et al.l {20021): ICharlton et alJ 
fcOOfl: IZonak et alJ j2004l) : IStocke et alJ <2004): 
iDahleml <2005l): ITumlinson et alJ 
U 120051) : iMasiero et al.1 <2005l) and 



Zonak^tal 
iBouche et alJ (120041 
J2005I): iKeenev et al 
I.Tenkins et al.l (|2005h . 

Papers supporting (b) incl ude (in ch ro nologi - 
cal order): iMorris et al.1 (Il99ll) : IMorris et all 1 199? ): 
Morris fc van den Bergbl (119941)7 iMo fc Morris! \l99i\: 



Wevmann et al.l (119951 ): iDinshaw et al.1 Jl995h:IStocke et all 
Jl995h: IShulletalJ <1996ft: iBahcall et al.l <1996ft : 
iRauch et alJ Jl99fJ): iLe Brun et al.1 <ll996l): iBowen et al 



^99); lvaj_Gc£kom_et al.l 
1 1998): |Weyman n et al] 



jDmshav^j^d 



Geller 



vandcn Berk et al 
Penton et alJ (l2000f 



J1999D: llmpev et al[ 

-OOd): iPenton et alJ <l2000l); iPenton et al 

fhl l i iDi) et al.M2002l):lTripp et all (l2002l) : jMcLin et al 

iMannind J2002D : iMcLin et alTj2002l) : Bowen et al 

iMa: 
1 i 



iRosenberg et al] J2003ft : iMannind <2003lk 
IPenton et alJ <2004|) : i Bregman et alJ <2004 : ISembach et alJ 
J2004h: ICote et all I2OO5I) : iDanforth fc Shull J2005h and 
IChen et al.l | |200SD . 



Finally to complete the reference list, some of the 
theoretical papers from the la st ten years w hi ch ar e 
relevant for this debate are: [Hernquist et al] il996 



Zhang et al.l (ll997l) : [Cen et al.1 (IT998t) :'lTheuns et al] ll998 
.Dave et alJ (Il999 | ) : ben fc Ostriked lll999l): iDaye fc Trippl 
J200ll): ISchave! J200l|): iMcDonald et alJ i200a)T 



IScannapieco et all J2OK ): 1viel et alJ J2002l): iTheuns et alJ 
J2002I): I McDonald ct al. J2002I): iNagamine et al.l j2004al);' 



iNagamine et al.l (l2004bF: iFuiita et alJ J2004I) : lAguirre et alJ 
J2005h: ITumlinson fc Fand (120051) : IScarmapiecol i2005l): 
iFurlanetto et al.l J2005I) : IViel et al.1 J2005h and IJena et alJ 



(200 

Regardless of whether either hypothesis is correct, we 
emphasize that even if one (naively) believes a clear answer 
can be obtained (either (a) or (b) above), this answer will 
be a function of the neutral Hydrogen column density of the 
absorbers. For Nhi> 10 21 cm -2 , one is probing a column 
density of material comparable to that of the disk of our 
galaxy, and hence all galaxies, to varying degrees, should be 
expected to contribute to the population of such absorbers. 
At column densities Nhi~ 10 12 cm -2 , one is getting close 
to the neutral hydrogen content of the expected fluctuations 
within voids, and we expect that a more heterogeneous set 
of causes could produce such absorption. 

The increasingly sophisticated models of galaxy forma- 
tion and evolution now available (see references above) sug- 
gest that if we could make a perfect census of the gaseous 
and luminous constituents of a large volume of the low red- 
shift Universe, we would find that both the gas (quasar ab- 
sorption line systems) and stars (galaxies) trace the same 
fundamental structures (whose course of formation and dis- 
tribution was set by the early perturbations of the distribu- 
tion of dark matter in the Universe). These same models, 
however, also indicate that a wide range of detailed rela- 
tionships between the gas and galaxies should be observed. 
The details of the structures will depend on an equally wide 
range of interesting astrophysics, including the process of 
star formation and "galaxy feedback" and we hope that the 
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Figure 1. Illustrative image of the two lines-of-sight (LOS) to 
the QSOs LBQS 0107-025A and B. The region shown is 7.1'x6.2', 
with North up and East to the left. The two QSOs discussed in 
this paper are label led, as is a third QSO (labelled C) to the North 
iYoung et al.l200ll) . A few of the measured galaxy redshifts in the 
field are labelled (drawn so the decimal point is approximately 
above the relevant galaxy) and the projected spatial offset and 
velocity difference to detected absorbers in either of the LOS are 
marked for two example galaxies. See the text for a discussion. 
The image shown was made from a 90 second exposure without a 
filter using the Palomar Observatory 5m telescope and COSMIC 
imager/multi-object spectrograph on September 22, 1995UT. 



approach of this paper can lead to a better understanding 
of this astrophysics. 

Observationally comparing the distribution of galaxies 
and quasar absorption line systems shows some of the com- 
plexity of the phenomena at work. This can be illustrated 
by considering Figure Q In this figure we show an image of 
the lines-of-sight (LOS) towards the QSOs LBQS 0107-025A 
and B. As discussed later in the paper, for this field we have 
samples of galaxies and Ly-a absorption line systems whose 
locations in this common volume of space can be compared. 
In the figure we have selectively labelled illustrative cases 
of 'associated' galaxies and absorption lines. For two of the 
galaxies labelled with a redshift, arrows indicate where there 
is a detected absorption line in both of the QSO spectra. For 
those absorber-galaxy pairs, the projected spatial separation 
and velocity difference between the galaxy and the absorber 
are also provided. 

We consider a few of these cases below: 

(i) A bright, multi-component object lies between the two 
QSOs at a redshift of 0.2272. There is absorbing gas seen at 
nearly this same redshift in both QSO LOS, and in one of 
them, CIV absorption is also detected. This then could po- 
tentially be interpreted as a straightforward case of a large 
gaseous halo around a bright galaxy. The complex morphol- 
ogy of the galaxy could indicate an ongoing merger. 

(ii) At a distance of 350 kpc to the south and east of the 
two QSO LOS is another bright morphologically-complex 
galaxy at a redshift of 0.1145. This galaxy is also close in 
velocity to absorption seen in both LOS. It is instructive 
though to compare the small visible extent of the stars in 
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the (large) galaxy with the distance to the two QSO LOS. 
Smooth, undisturbed gaseous halos of this size seem im- 
probable, although it is one of the goals of this paper to 
investigate this statistically. 

(iii) Generally to the north of the two QSOs are a collec- 
tion of 5 galaxies marked with redshifts near 0.20. All 5 of 
these galaxies are within 500 km s" 1 (and also within a pro- 
jected distance of 500 kpc) of an absorber seen in the LOS 
to LBQS 0107-025B. No absorption is detected in the other 
QSO LOS. This collection illustrates the common difficulty 
of picking out one individual galaxy to tie to one absorber. 
Although there is indeed a single 'closest' galaxy, it seems 
rather arbitrary to claim that it is the sole cause of any 
absorption. 

(iv) Finally, a collection of six galaxies, all within 500 km 
s _1 of each other at a redshift around 0.24 are marked sur- 
rounding the LOS to LBQS 0107-025A. Despite their prox- 
imity to the LOS, and the case above where a group of galax- 
ies can be associated with an absorber, no absorption is seen 
at this redshift in either of the two QSO LOS. While one 
could speculate that low column density absorbers might be 
fragile entities that are destroyed in a dense environment, 
there is no good independent evidence at present that this 
galaxy 'group' has crossed such a threshold. 

These examples demonstrate the great variety of ap- 
pearances the underlying relationship between the absorbers 
and galaxies presents, along even one LOS. This richness 
should be kept in mind when interpreting the results de- 
rived from the ensemble of data we consider in this paper. 

The outline of the paper is as follows: in §[5] we describe 
the various data sets used in this paper. In § [3] we analyse 
the data, and in §0]we present our conclusions. 



2 DATA 

Our observational program would ideally be designed to al- 
low a statistical comparison between the absorption gas de- 
tected in Hubble Space Telescope (HST) UV spectroscopy 
between redshifts of and 1 and the galaxies within a 
cylindrical region at least 3 Mpc 1 in radial extent away 
from the QSO Line of Sight (LOS). In practice, what we 
achieved was an absorption line s ample taken from th e HST 
Quasar Absorption Line Survey (|Bahcall et alJll99.1 Il99ft 
Ijannuzi et al.lll998t) . an HST Key Project during cycles 1-4, 
with varying redshift coverage, and galaxy samples from a 
variety of observatories. The main one reported here was 
obtained using the Multi-Object Spectrograph (MOS) at 
the Canada- France-Hawaii Telescope (CFHT). Some addi- 
tional galaxy redshifts have been obtained using the COS- 
MIC spectrograph at Palomar, and gleaned from the lit- 
erature. We are in the process of supplementing this data 
set with additional observations using other facilities. Al- 
though we do not claim that any of these galaxy samples 
are 'complete', we have endeavoured to ensure that they 
all come with known selection functions, allowing a statisti- 
cal analysis to be performed. Completeness in this context 

1 Throughout we assume Hq=70 km s _1 Mpc -1 , f2=0.3 and 
A=0.7. Impact parameters are reported in physical rather than 
comoving coordinates. 



would be rather an illusory concept, as there will always be 
weaker lines and less luminous or lower surface brightness 
galaxies that are missed by any su rvey. This point has been 
made by many author s including Linder (1998, 2000) and 
iRosenberg et alJ ll2003lh 

2.1 Sample definition 

The fields observed with M OS on the CFHT wer e chosen to 
include Key Project QSOs lljannuzi et aljfl998f) with QSO 
redshift greater than 0.5 and Faint Object Spectrograph 
(FOS) higher dispersion spectroscopy covering from redshift 
0.34 to 1 (i.e., FOS gratings G190H and G270H). While UV 
spectroscopy down to z=0 is obviously desirable (and has 
been obtained with the HST along a limited number of sight 
lines) the KP LOS including coverage to z=0 (observations 
made with the FOS G130H grating) are of QSOs with red- 
shifts too low for efficient CFHT MOS followup. In Table □ 
we list the QSOs selected. The columns include the QSO 
RA, Dec, redshift and V band magnitude, followed by the 
date of the HST UV spectroscopy for the 3 FOS gratings 
used in this paper, followed by the number of CFHT MOS 
masks observed (see § 12.41 and the date of the MOS obser- 
vation. 

There is one additional field observed as part of the 
CFHT MOS sample which was not part of the HST Key 
Project. HST UV spectra of the QSO pair LBQS 0107-025A 
and LBQS 0107 - 025B were studied bvlDinshaw _et ail <199Eft , 
iDinshaw etafl dl997fl and I Young et alJ <200ll) . This LOS 
was included in the CFHT MOS target list because of the 
presence of two QSOs bright enough for reasonable S /N HST 
UV spectroscopy within the MOS field of view. An analysis 
taking full advantage of this will be performed in a later 
paper, but for the present paper we generally use them as 
two 'separate' LOS. We have checked whether removing one 
of the LBQS 0107-025 QSOs from the sample makes any 
significant difference to our numerical results or conclusions, 
and have found that it does not. 



2.2 HST Quasar Absorption Line Sample 

The complete list of selectio n criteria for the KP quasars 
is discussed in section 2.1 of Ijannuzi et alJ lll99ct) . We re- 
view here only the criteria used by the KP that might pro- 
duce a bias for our study. First is their practical decision 
to select targets with bright V-band magnitudes and well 
determined redshifts in order for their targets to have pre- 
dictable UV luminosities that, in the absence of high column 
density absorption systems along the LOS, would enable a 
high signal-to-noise ratio UV spectrum to be obtained in a 
modest exposure time with HST. A sample biased toward 
objects bright in V-band apparent magnitude might be more 
inclined to include lensed quasars, i.e. objects with a fore- 
ground over-density along the LOS. In practice, the redshifts 
of the quasars included in the KP are low, and therefore 
less likely to be lensed. Second, a subset of the quasars in- 
cluded in the KP sample were not randomly selected, but 
were chosen for a GTO program of J. Bahcall because they 
were behind foreground objects whose absorption properties 
were to be investigated. These include AGN like Markar- 
ian 205 (with the galaxy NGC 4319 very close to the LOS 
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Table 1. CFHT MOS Sample Definition 



OBJECT 


R.A.(J2000) 


Decl.(J2000) 


z 


V 


G190H 


G270H 


G160L 


Masks 


Date 


NAB 0024+22... 


00 


27 


15.4 


+22 41 59 


1.118 


16 


6 


07/22/92 


07/22/92 


05/11/91 


1 


Jul-95 


LBQS 0107-025A 


01 


10 


13.1 


-02 19 52 


0.960 


17 


9 


02/05/97 


12/26/96 




1 


Jul-95 


LBQS 0107-025B 


01 


10 


16.2 


-02 18 50 


0.956 


17 


3 


02/12/94 


12/26/96 




1 


Jul-95 


PKS 0122-00... 


01 


25 


28.9 


-00 05 57 


1.07 


16 


7 


11/09/91 


09/11/91 




2 


Nov- 9 7 


3C 57... 


02 


01 


57.2 


-11 32 34 


0.67 


16 


1 


11/01/95 


11/01/95 




1 


Nov-97 


3C 110... 


04 


17 


16.7 


-05 53 46 


0.773 


15 


9 


03/16/95 


03/16/95 




3 


Nov-97 


NGC 2841 UB3... 


09 


19 


57.7 


+51 06 10 


0.553 


16 


5 


03/04/92 


03/04/92 


03/04/92 


2 


Nov-97 


0959+68W1... 


10 


03 


06.8 


+68 13 18 


0.773 


15 


9 


11/12/92 


11/12/92 


11/12/92 


1 


Nov-97 


4C 41.21... 


10 


10 


27.5 


+41 32 39 


0.613 


16 


9 


09/12/92 


10/12/92 


10/12/92 


2 


Nov-97 


3C 334.0... 


16 


20 


21.8 


+17 36 24 


0.555 


16 


1 


03/04/92 


03/04/92 


03/04/92 


2 


Jul-95 


PG 1718+481... 


17 


19 


38.3 


+48 04 12 


1.084 


11 


7 


05/13/93 


05/13/93 




1 


Jul-95 


PKS 2145+06... 


21 


48 


05.5 


+06 57 39 


0.99 


16 


5 


10/22/91 


10/22/91 




3 


Jul-95 


3C 454.3... 


22 


53 


57.8 


+16 08 53 


0.859 


16 


1 


11/15/91 


11/15/91 


11/15/91 


2 


Jul-95 


PG 2302+029... 


23 


04 


45.0 


+03 11 46 


1.052 


15 


8 


05/10/94 


02/08/94 




2 


Nov-97 


PKS 2340-036... 


23 


12 


56.6 


-03 22 27 


0.896 


16 





10/28/91 


01/18/93 


01/18/93 


2 


Nov-97 


PKS 2344+09... 


23 


16 


36.9 


+09 30 46 


0.677 


16 





07/10/92 


07/10/92 


07/10/92 


2 


Nov-97 



to Mkn 205, prod ucing absorption i n the spectrum of th e 
background AGN; IWeedmanl (119701) : IBahcall et all (I1992TH 
and 3C 232 (with the galaxy NGC 3067 in the foreground) . 
Of the objects we selected from the KP catalogues for obser- 
vation with CFHT, only NAB 0024+22 was not a randomly 
selected target by the KP. This quasar was part of the GTO 
program of Bahcall and was selected because it was discov- 
ered as par t of a search for qua sars with LOS near a cluster 
of galaxies iBahcall et al.llT973l) . Abell 31, with a redshift of 
0.160, is along the LOS to this quasar. The cluster redshift 
is smaller than the range we investigate along this LOS. 

We had two additional criteria for selecting our sample 
for the CFHT observations. The fields needed to be observ- 
able during our scheduled CFHT observing runs. We chose 
to focus on LOS with quasars with redshifts greater than 
0.5, in order to maximize the efficiency of our redshift galaxy 
survey. 

With the previously mentioned exception of 
LBQS 0107— 025A and B, the sample of absor ption lines we 
use is d rawn from the KP ca talogu e papers. [Bahcall^^tlJ 
lll993ft . IBahcall et all lll996ft . and Ijannuzi et alJ <199fl) . 
The KP Ly-a absorption line catalogue was constructed to 
include all features with a signifi cance level grea t er tha n 
4.5a (as defined and discussed bv ISchneider et all (^993)). 
This corresponds to a variable equivalent width limit along 
each LOS. Plots showing these equivalent wi dth limits as a 
functi on of wavelength are given in Figure 2 of lJannuzi et alJ 
(1998). With the modest spectral resolution of the FOS 
(R=1300; approximately 230 km s _1 ), the Ly-a absorption 
features were generally unresolved and conversion from the 
observed equivalent width to column density in neutral 
hydrogen must be done with some assumption about the 
Doppler-b parameter for the lines. For the subset of lines 
that were broader than the instrumental profile, their width 
probably indicates the blending of multiple absorption 
features, and not a line profile that can be fit to deduce the 
column density of the absorber. We note that approximately 
10% of the detected Key Proje ct Ly-q absorbers we re in 
fact resolved, and, as noted by Ijannuzi et alJ dl99Sf) . are 
a sign of possible clustering of absorbers on these velocity 



scales. This results in an additional contribution to the 
uncertainty in the derived column densities. 

For LBQS 0107- 025A and B we h ave chosen to 
use the line lists from lYoung et alJ ll200ll) , supplemented 
with l ines from the GHRS observations of iDinshaw et alJ 
(Il997l) . (I n detail this m e ant a dopting the lines from Ta- 
ble 5 of IDinshaw et all Jl99Tt> . and Tables 3 and 5 of 
lYoung et alJ 12001 ) for L BQS 0107-025A, and from Table 5 
of TDinshawel^n' T ll997t) . and Tables 4 and 6 of lYoung etafl 
feOOll) for LBQS 0iq7-025 B.) We shifted the wavelengths 
of the | D inshaw et alJ (Il997f) by -1.72 A as recommended by 
lYoung et alJfcOOll) . 

We want the subset of absorbers caused by interven- 
ing, rather than assoc iated, absorption systems. We follow 
IWevmann et alJ dl998l) and exclude any line from our sample 
with a difference of less than 3,000 km s _1 from the systemic 
velocity of the quasar. Note that while this is a traditional 
choice of velocity range for excluding associated systems, it 
will not exclude all associated systems. In particular, one 
of our quasars, PG 2302+029, includes at least one asso- 
ciat ed system wit h an e j ection veloci t y of o ver 50,000 km 
s -1 J Jannuzi et alJ il99r3) : ISabra et alJ ibooaft ; Jannuzi et al. 
2006, in preparation). We have excluded the broad system 
at z — 0.695 from the sample used in this paper. 

All of the line lists we use have complete identifications 
with the exception of the spectrum of PG1718+481. The 
effects of the incomplete identification of the PG1718+481 
absorbers are discussed in § 13.41 where it is shown they are 
negligable. 



2.2.1 Comments on Individual Lines of Sight 

Detailed notes about the HST FOS spectra along the indi- 
vidual LOS we are studying, and the resulting absorption 
line lists, are inclu ded in the KP catalogue papers and in 
lYoung et alJ fcOOllU for the LBQS 0107-025 LOS). Here we 
comment on the special properties of a few of the fields. 

As previously mentioned, NAB 0024+22 was discovered 
durin g a search for quas ars near or behind clusters of galax- 
ies (IBahcall et alJll973f> . Abell 31, with a redshift of 0.160, 
is foreground to the redshift range covered by the FOS spec- 



The association between gas and galaxies 5 



troscopy (0.2950 to 1.0998). Intervening absorption line sys- 
tems including absorption by CIV are at z = 0.4069, 0.4830, 
0.8196, and 1.1102. The SNR of the KP spectrum of this 
quasar is not as high as the typically observed object, re- 
sulting in a smaller sample of detected absorbers given the 
rela tively long redshift p ath ob served (Sz = 0.8). 

iBergeron fc Boissel lll99ll) reported an emission line 
galaxy at z=0.791 in the LOS to PKS 2145+06, which was 
later identified as being associated with an extensive metal 
line systems in the F OS UV spectrum obtained by the KP 
iBergeron et alll994h . This z = 0.791 system is high enough 
column density to produce a Lyman limit system starting 
at 1633A , preventing the observation of Ly-a absorption 
along this LOS below z~0.3. 

2.3 Imaging 

To compare the distribution of the galaxies along the same 
LOS that the FOS spectroscopy maps the gas distribution 
requires knowledge of the positions and redshifts of galax- 
ies in the fields of our quasar sample. Images of the quasar 
fields can provide the positions and magnitudes of the galax- 
ies. The spectroscopy that allows the determination of the 
galaxy redshifts requires the imaging catalogue as input for 
the design of multi-slit masks. 

For the CFHT MOS, the images used to generate the 
multi-slit mask designs are obtained during the same ob- 
serving block as the spectroscopy (i.e., the instrument has 
not been removed from the telescope between the imag- 
ing and spectroscopy observations). Therefore we obtained 
images for our fields during our two observing runs. The 
images were reduced shortly after they were obtained and 
catalogues of galaxies generated to guide the design of the 
multi-slit masks. After the observing runs were completed, 
the images were re-reduced, the photometry and astrome- 
try of the images calibrated, and final catalogues were con- 
structed. Since the main purpose of the scheduled observing 
time was to perform the spectroscopy, only R— band images 
deep enough to detect the galaxies we wished to target for 
spectroscopy were obtained. It would be desirable to obtain 
images including additional band-passes and with enough 
photometric standard observations to improve the accuracy 
of our photometry. 

In the following subsections we describe the imaging 
observations, the astrometric and photometric calibration of 
the images, and the measurement of the galaxy positions and 
J?— band magnitudes in the fields of our targeted quasars. 

2.3.1 Observing 

We had two observing runs with CFHT MOS (July 29 - 
August 1, 1995 UT and November 29 - December 3, 1997 
UT). The date each field was imaged is listed in Tabled The 
instrument and telescope properties were significantly differ- 
ent during the two runs and we will describe the particulars 
by run. 

During the 1995 observing run, the 29th of July UT 
was photometric and 600 second R— band exposures were 
obtained of six quasar fields. Two and 10 second exposures 
of the photometric sta ndard stars in the field of NGC 7790 
lOdewahn et alJ lll992l) were also obtained to allow a pho- 



tometric zero-point to be determined for the images. Un- 
fortunately, the proper baffles were left off of the telescope 
when MOS was installed on the telescope prior to our run. 
As a result, bright arcs of light, produced by stars outside 
of the FOV, cut across the images and effectively increased 
the sky level over approximately one third of the FOV of the 
images. The amount and distribution of the scattered light 
was, naturally, pointing dependent. The mixture of sharp 
and diffuse features in the scattered light made it difficult 
to completely remove during the processing of the images. 
This complicated the determination of the sky background 
around objects, and the subsequent generation of the galaxy 
catalogues necessary for the construction of the multi-slit 
masks. The two most strongly affected images were those of 
LBQS 0107-025AB and PC 1718+481. To increase the area 
that could be mapped for galaxies, we obtained two addi- 
tional images of each field, each offset in position from the 
others, and combined the images of each field after masking 
the regions with the worst scattered light in each of the expo- 
sures. The combined images allowed us to produce improved 
measurements of the galaxies in these fields, but there are 
still regions in each image where we would not be able to 
detect galaxies and objects whose measured properties were 
strongly affected by scattered light. The imaging on 30 and 
31 July 1995 UT was not done under photometric condi- 
tions, so the flux zero-point for the combined image was set 
to match the imaging on the 29th. The CCD in MOS during 
this run was the Loral-3 with pixels of approximately 0.313" 
in size and a gain of 1.9 electrons per ADU. 

During the 1997 observing run, the 29th of Novem- 
ber UT was photometric and 600 second exposures of four 
quasar fields were obtained. Two and 10 second exposures 
of the photometric standard stars in the field of NGC 7790 
were also obtained. The correct baffles were on the tele- 
scope. The nights of UT 26 Nov, 30 Nov, and 1 Dec, during 
which 600 second exposures of 5 additional quasar fields 
were obtained, were not photometric. For two of the fields 
taken during these nights, two 600 exposures were obtained 
to compensate. The CCD in MOS during this run was the 
STIS-2 with pixels of approximately 0.44" in size and a gain 
of 4.52 electrons per ADU. 



2.3.2 Reductions 

The images were reduced and analysed twice. During the 
observing run a basic reduction was performe d with I RAF 2 
(bias subtraction and flat fielding). FOCAS dValdeslll989h 
was used to generate a catalogue that could be used for the 
identification of galaxies and the generation of the multi- 
slit masks. Following the observing runs, the images were 
processed in a more relaxed manner, astrometric and pho- 
tometric calibrations were determined, and measurements 
of the object positions and magnitudes completed. In the 

2 IRAF, the Image Reduction and Analysis Facility, is a general 
purpose software system for the reduction and analysis of astro- 
nomical data. It was written and supported by programmers in 
the Data Products Program of the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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Table 2. Image Quality 



Object 


UT Date 


DIQ 

u 


Comment 


NAR 0024+22 
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1. 
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T ROS 0107 09 C »AR 


9Q Tul Q^i 
j til yj 


1. 


10 








Sf) Tul Q^i 


0. 


80 


2 


exposures 


PK^ 01 99 on 


9Q Near Q7 


0. 


99 






3C 57... 


30 Nov 97 


0. 


.85 






3C 110... 


29 Nov 97 


1. 


22 






NGC 2841 UB3... 


26 Nov 97 


1. 


23 


2 


exposures 


0959+68W1... 


29 Nov 97 


0. 


95 






4C 41.21... 


1 Dec 97 


1. 


14 






3C 334.0... 


29 Jul 95 


1. 


15 






PG 1718+481... 


29 Jul 95 


0. 


68 








31 Jul 95 


0. 


59 


2 


exposures 


PKS 2145+06... 


29 Jul 95 


0. 


66 






3C 454.3... 


29 Jul 95 


0. 


.88 






PG 2302+029... 


29 Nov 97 


0. 


80 






PKS 2340-036... 


30 Nov 97 


1. 


17 






PKS 2344+09... 


26 Nov 97 


1. 


01 


2 


exposures 



following subsections we describe this second round of re- 
ductions, which yielded the measurements presented in this 
paper. 

The final processing of all of the images was performed 
using cedproc and other tasks distributed with v2.12 of 
IRAF. Bias subtraction and flat fielding were the main pro- 
cessing steps. While only single exposures were obtained for 
most fields, effort was made to identify saturated pixels and 
defects caused by cosmic rays, scattered light, satellite trails, 
etc., and to map the locations of affected pixels in a bad 
pixel mask for each image. These masks were used during 
the generation of the catalogues, described below, to flag 
objects whose properties could not be well measured. 

2.3.3 Astrometry 

In order to generate relatively accurate positions for the 
galaxies whose properties are presented in the electronic Ta- 
bles 3-18, we needed to determine for each of the images 
an accurate world coordinate system. This was generated 
using software in IRAF (tasks ecse twes and mscc match in 
the mscred package; as described bv lValdesl (12002ft 1 and the 
USNO A2. catalogue of the positions of objects in each of 
the fields iMonet et al.l fl998h . A limited magnitude range 
was used of USNOA2.0 reference stars in order to avoid the 
documented relative offset in astrometry for this catalog as 
a function of magnitude. Using the IRAF task msccmatch 
(in v4.7 of the mscred package), a fourth order mapping be- 
tween the images in x,y pixel coordinates to RA and DEC 
of the reference stars was determined. The RMS residuals of 
the fit were 0.3 to 0.4" for each of the images. At the time 
these solutions were determined, the USNO A2.0 was the 
best catalogue available. Although catalogues with improved 
positio ns are now available (e.g. USNO-B1.0, IMonet et alJ 
J2003ft or GSC2.2 3 ), a trial new solution using the USNO- 

3 The Guide Star Catalogue-II is a joint project of the Space 
Telescope Science Institute and the Osservatorio Astronomico 
di Torino. Space Telescope Science Institute is operated by the 



B1.0 did not result in a significantly improved solution, so 
we have retained our original fits. The residuals to our so- 
lutions were largest around the border of the images and 
we suspect that the uncertainties in the positions of objects 
around the edge of each field should probably include an 
additional 0.1 to 0.2" of systematic error. Our measured po- 
sitions of objects given in this paper are in the ICRF and 
J2000. 



2.3.4 Photometry and Catalogue Generation 

We did not obtain multiple standard star observations nor 
observe in a range of colours that would allow us to de- 
termine our own extinction solutions and colour terms. This 
fundamentally limits the accuracy we can obtain in our pho- 
tometry. However, our needs in the current study only re- 
quire photometry accurate to a few tenths of a magnitude. 
Using a representative Mauna Kea extinction curve deter- 
mined from photometric observations during other runs and 
our observations of the stars in NGC 7790, which h ave been 
calibrated in the R— band bv lOdewahn et all il992t) . we were 
able to determine zero points for the 10 fields imaged during 
our photometric nights. 

There were five fields observed only on non-photometric 
nights. We have determined an approximate zero-point for 
these images as follows. First, for the 10 fields we calibrated 
with our measurements of NGC 7790, we compared our mea- 
sured magnitudes to the red magnitudes of the GSC2.2 for 
objects in common. The pass-bands, of course, are not the 
same, and we would not expect the magnitudes to be iden- 
tical, but we were gratified that all of the fields observed 
under what we believe to have been photometric conditions 
yielded similar offsets (mean offset of 0.1 of a magnitude) 
and RMS residuals (0.15 to 0.25). We then assigned to each 
of the five fields (all from the 1997 run) observed on non- 
photometric nights a zero-point based on the calibration de- 
termined from the photometric night of that run. We then 
compared these five fields to the GSC2.0. Two of the fields, 
3C 57 and PKS 2340—036, had mean differences in the same 
range as those fields observed under photometric conditions. 
Three of the fields; NGC 2841 UB3, 4C 41.21, and PKS 
2344+09; had significantly larger initial offsets (0.43, 0.47, 
and 0.62 magnitudes). The RMS residuals, however, were 
similar to those of all the other fields. We therefore adjusted 
the zero-points for these fields so that their mean offset when 
compared to the GSC2.0 stars would also be 0.1 of a mag- 
nitude. This adjustment is consistent with what is evident 
from inspecting the images (and number counts), namely 
that these exposures taken with some clouds did not reach 
the same depth (in the same exposure time) as our other ex- 
posures. A systematic uncertainty of 0.2 mags should prob- 
ably be assigned to our reported magnitudes for observed 



Association of Universities for Research in Astronomy, for the 
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NAS5-26555. The participation of the Osservatorio Astronomico 
di Torino is supported by the Italian Council for Research in 
Astronomy. Additional support is provided by European South- 
ern Observatory, Space Telescope European Coordinating Facil- 
ity, the International GEMINI project and the European Space 
Agency Astrophysics Division. 
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galaxies, but we feel we have successfully placed our obser- 
vations on a common relative photometric scale. 

Using our magnitude zero-points for each image and 
measurements of the variance in the sky brightness of each 
image, we have confirmed that the five sigma detection 
thresholds for point sources in all of our images is a magni- 
tude of 24.9 or greater, well below our selection of targets 
for spectroscopic follow-up. 

For each field th e object catalogues were generated us- 
ing SExtractor 2.2.2 jBertin fc Arnoutsl[l993l run with the 
minimum detection area, Gaussian convolution filter, and 
signal above sky threshold optimized to detect and measure 
all of the objects for which spectroscopic data was obtained. 
We measured eight fixed aperture magnitudes (starting at 
2" and increasing to 10" in diameter) and the SExtrac- 
tor MAG_BEST, which is similar to the Kron total magni- 
tude llKrorJl980l) . Results for our spectroscopically observed 
galaxies are included in the electronic Tables 3-18. 

As previously discussed, for some of our fields scattered 
light affected our images. Even without these problems, any 
catalog will suffer increasing incompleteness for galaxies of 
fainter total magnitude and/or surface brightness. As a test 
of the completeness of our imaging catalogs, artificial galax- 
ies were inserted into each image (using the IRAF artdata 
task) and SExtractor used to attempt to recover the objects. 
While not a perfect measure of the robustness of our cat- 
alogs, the 50% completeness limit of our catalogs for each 
quasar field is typically 24.4 (R— band) and greater than 24.0 
for all fields. The fields most seriously affected by the arcs of 
scattered light could be missing even bright galaxies at the 
5% level, but visual inspection indicates that the fraction is 
likely to be much smaller. 



2.3.5 Bounded Restframe B-band Magnitude 

We would like to have well measured luminosities for each 
of the galaxies expressed in a common rest-frame band. Un- 
fortunately, we do not have well measured colours or spec- 
trophotometry that would allow an accurate determination 
of the over-all spectral energy distributions (SEDs) of the 
galaxies in our sample. This prevents the assignment of the 
proper A'— correction for each galaxy in our sample. 

We can, however, constrain the range of rest frame 
B— band luminosities possible for each galaxy using our 
R— band measurements, the redshifts determined from 
our spectroscop y, and the galaxy S EP templates and 
A'— corrections of lColeman et all ([j.980). For each galaxy we 
determined the range of possible K~ corrections and intrin- 
sic B — R colours and tabulated the resultant maximum and 
minimum rest-frame B— band magnitudes. These are listed 
for each galaxy in the electronic Tables 3-18. We assumed 
no evolution or intrinsic extinction. While more involved 
procedures might, with other data sets, provide probabil- 
ity distributions for the rest B— band luminosities of each 
galaxy, the procedure described above allows us to make the 
gross comparisons between sub-samples that are reasonable 
given the other uncertainties in our measurements. In other 
words, those galaxies that might be the intrinsically more 
luminous are separable from those that could not. The ob- 
served redshifts and R— band magnitudes are the dominant 
constraints. 



2.4 MOS Spectroscopy 

2.4-1 Observing 

The majority of the galaxy redshifts described in this paper 
were obt ained using the CFH T Multi-Object Spectrograph 
(MOS), l)Le Fevre et alll994T) . during observing runs in July 
1995 and November 1997. As described above, imaging cat- 
alogs were prepared in real time, from which masks were 
designed. For July 1995 the Loral 3 CCD was used along 
with the O300 grism, yielding roug hly 3. 5 A per pixel. No 
wavelength blocking filter was used, yielding a wavelength 
coverage that varied depending on the location of the slit 
within the field of view, but reasonable signal-to-noise (S /N) 
was generally obtained for wavelengths between 5000 A and 
9000 A. For the November 1997 run, the STIS 2 CCD was 
used with the same grating, yielding a dispersion of 5.2A per 
pixel, for both runs, masks were designed with a 1.5" slit 
width, corresponding to 4.8 pixels in July 1995 and 3.4 pix- 
els in November 1997. Some of the observations were taken 
through cloud. A flux standard star was observed in order 
to allow an approximate flux calibration of the data and to 
allow us to remove the instrumental signature, but because 
of the above mentioned cloud, the fluxes recorded could be 
low. 

Spectroscopic multi-slit masks were designed using the 
galaxy catalogs generated during the observing run. Slits 
were assigned to galaxies with a magnitude ranking. Slits 
had a default length of 10" and width of 1.5". After this, the 
assignment algorithm checked whether it was possible to add 
further objects by sliding slits sideways by a limited amount. 
Finally all slits were expanded in length to fill the array and 
allow better sky subtraction for some objects. Depending 
on the galaxy surface density and distribution in the field 
between 40 and 47 galaxies had slits assigned to them per 
mask in this manner. Between one and three masks per LOS 
were taken. 

Because of the varying conditions mentioned above, the 
redshift success rate for objects which were assigned slits 
varies from mask to mask. The average success rate was 0.5, 
with a nadir of 0.07 for a mask taken through thick cloud, 
and a pinnacle of 0.89 for one particularly good mask. 

It was never our intent to obtain a sample complete to 
any given magnitude, but merely to obtain a statistically 
useful sample of the galaxies around the LOS. By compar- 
ing with the photometric catalog, we can reconstruct our 
completeness for each field as a function of magnitude, as 
will be discussed below. 

2.4-2 Reduction 

The spectroscopic data reduction using IRAF followed the 
standard procedure for MOS spectrosco py used for example 
in the CNOC2 project <Yee et al.lE'o OO'). The data were bias 
subtracted and trimmed. If necessary, bad columns were in- 
terpolated over. Repeat exposures of the same mask were 
combined and cosmic rays removed. Because the spectral 
extraction involves summing over several columns, flat field- 
ing to correct for pixel-to-pixel variations in CCD sensitivity 
was not necessary (and would have added noise) . Individual 
apertures were then summed to generate ID spectra, simul- 
taneously subtracting the sky using adjacent regions along 
the slit. These spectra were wavelength calibrated using arc 
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spectra obtained at the same position on the sky (to min- 
imize the effects of flexure), and the spectra were approxi- 
mately flux calibrated using observations of a standard star 
taken through the same slit width. 

2-4-3 Redshift Determination 

The driving goal for the redshift measurements was to ob- 
tain redshifts that could all be placed in the same reference 
frame, in our case heliocentric. To do this, all the galaxy 
redshifts determined in different ways were tied together in 
the same reference frame, and then this frame was tied to 
that of the absorption line lists, which were adjusted to the 
heliocentric velocity frame as described by Jannuzi et al. 
1998. 

All galaxies were cross correlated with two absorption 
line templates using the IRAF XCSAO task. One template 
was the 'fabtemp97' one which comes with the XCSAO pack- 
age, and the other was a high-quality wide- wavelength cov- 
erage spectrum of the early-type galaxy NGC 4889, kindly 
provided by B. Oke (private communication, template avail- 
able on request from S. Morris). As usual, the cross corre- 
lation process involved continuum fitting and also filtering 
of high and low Fourier components. Redshifts for emis- 
sion line objects were obtained more interactively using the 
IRAF RVIDLINES task inside the RV package. The results 
of these routines were then checked by eye, and each mea- 
surement flagged as successful or failed. This somewhat sub- 
jective procedure was necessary in order to deal with the 
combined effects of zero order contamination, poor sky sub- 
traction, residual cosmic rays and other reduction problems. 
In practice, for the cross correlation measurements, there 
was a fairly clear threshold in R-value reported by the rou- 
tine above which between measurement were judged success- 
ful. Despite all of the above, it is still of course possible for 
some objects in our redshift catalog to have 'catastrophic' 
redshift errors (i.e. redshifts based on spurious features of 
mis-identified features). We believe that the above proce- 
dure should have minimized the number of such errors, and 
would be surprised if there were more than a handful these in 
the redshift catalog, but it should be noted that the redshift 
errors do not take account of this possibility. 

Both the above routines automatically convert redshifts 
to a heliocentric frame. Nevertheless, as is well known, this 
does not mean that there are no systematic shifts in veloc- 
ity introduced by the different processes, and so objects with 
reliable redshifts determined from 2 or more of the above ap- 
proaches were used to measure and remove any such shifts. 
(Shifts of -49 km s^ 1 and +72 km s" 1 were measured for 
the NGC 4889 and emission line redshifts relative to the 
fabtemp97 measurements) . 

At this stage we are left with 641 galaxy redshifts in 
our catalog with their associated redshift errors. The success 
rate of obtaining redshifts from a given number of slits on a 
mask varied considerably depending on conditions and the 
magnitude range of objects on the masks, but for typical 
masks and conditions was around 60-70%. 

As described above, the photometry performed on the 
mountain for mask design was repeated after the observing 
run. The resulting photometric catalog was then re-matched 
with the redshift catalog, resulting in our dropping two of 
our measured objects because of large uncertainties in their 



positions and magnitudes. The CFHT MOS data therefore 
provided 639 new redshifts of galaxies in the targeted LOS. 

The objects with identified redshifts are listed in a set 
of Tables (3-18) available in the electronic version of the 
paper. The columns in the associated electronic table are as 
follows: 

(i) RA: The J2000 Right ascension of the galaxy. The 
relative positions for galaxies in a single quasar field should 
be good to a few tenths of an arcsecond. 

(ii) DEC: The J2000 Declination of the galaxy. The rela- 
tive positions for galaxies in a single quasar field should be 
good to a few tenths of an arcsecond. 

(iii) z: The redshift determined from our MOS spectrum 
of the galaxy. 

(iv) z crr : this is the error reported by the cross correlation 
routine, and as such should be considered a lower limit on 
the redshift error. 

(v) R-mag: R-band apparent magnitude (MAG_BEST 
from Sextractor). As discussed above, these are only accu- 
rate to ±0.2 magnitudes. 

(vi) B max : estimate of the maximum rest frame B band 
luminosity as discussed in § 12.3.51 

(vii) B m i n : estimate of the maximum rest frame B band 
luminosity as discussed in S 12.3.51 

2.5 Additional Galaxy Redshifts 

We also searched the literature to find additional galaxies 
with measured redshifts around the sample QSO LOS. As 
the typical redshift range for the SDSS and 2dF galaxy red- 
shift samples are not well matched to absorption line data 
used in this paper, we have not chosen to add in galaxies 
from those surveys at this stage. 

Our main method for finding additional redshifts was 
through the NED database 4 . We list the results of this search 
below. Since the analysis in this paper is statistical in nature, 
we have chosen only to add redshifts from other authors to 
our sample when they make a significant difference to the 
number of redshifts known for a given LOS. 

2.5.1 PKS 0122-00 

Galax ies in this LOS were ob served by both iBowen et alJ 
Jl997h . and IChen et alJ (1200 if) , but as only two additional 
galaxy redshifts would have been added to the sample, for 
simplicity these were not included. For t he three galaxies in 
common between our sample and that in lChen et alJ (l200lT) 
a mean velocity difference of 30 km s _1 is measured, with 
an RMS of 450 km s" 1 . 

2.5.2 PKS 2145+06 

A galaxy that seems to be associated with a Mgll absorber 
at z~0.79 i n the LOS has been s tudied in a number of 
papers, e.g. dChurchill et alJl2000allrl . This galaxy has not 
been added to the sample. 

4 This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 
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Table 3. Sample Table - CFHT MOS Sample Galaxy Properties - NAB 0024+22 



R.A.(J2000) 


Decl.(J2000) 


z 


z crr 


R-mag 


Bmax 


Bmin 


0:26:55.8 


22:41:02.3 


0.1587 


0.0003 


18.1 


-20.4 


-19.6 


0:26:56.7 


22:40:09.3 


0.3100 


0.0006 


20.2 


-19.9 


-19.3 


0:26:57.8 


22:41:30.6 


0.3469 


0.0007 


21.0 


-19.3 


-18.9 



2.5.3 PG 2302+029 

This LOS was surveyed by iBowen et alJ (120011b obtaining 
spectra for 24 objects. In this same field, we measured red- 
shifts for 42 objects. We had 6 galaxies in common. (For 
which we measure a reassuringly small mean velocity shift 
of -65 km s _1 with RMS 227 km s _1 , in line with our error 
estimates). This then gives us 18 additional redshifts which 
we add to our sample. These objects are included in (the 
electronic) Table 16. On a less positive note, it is interesting 
that 2 groups designing masks for the same LOS could end 
up with so little overlap in their samples. 

2.5.4 Q0107- 025A,B 

The region around the LOS to the QSO pair 
Q0107-025A,B was also observed with the COS- 
MIC spectrograph on the Palomar 5m telescope 
( http: / / www.astro.caltech.edu/palomar/200inch /cosmic/l l 
by Weymann et al. (private communication) . The measured 
redshifts for one of the masks observed was kindly made 
available to us by M. Rauche. This yielded 28 additional 
objects with redshifts. In this case there was no overlap by 
design, and so we are assuming that any systematic shifts 
between the two redshift sets is small. These objects are 
included in (the electronic) Table 4. 
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Figure 2. Pie diagram showing absorbers and galaxies in the 
LOS of NAB 0024+22. The top and middle panels show projec- 
tions of the galaxy distribution in RA and Dec, while the bottom 
panel shows the galaxy distribution in impact parameter as a 
function of redshift. Open circles represent galaxies, open trian- 
gles are absorbers detected in Ly-cr, filled triangles are absorbers 
seen in CIV. A star (if plotted) shows the location of the QSO. 
In all three panels, the vertical dashed line at z=0.34 shows the 
approximate location of the lowest redshift Ly-a detectable with 
the HST FOS G190H grating. In the bottom panel, the horizontal 
dashed line shows an impact parameter of 280 kpc in our adopted 
cosmology (see § 13.21 . the curved solid line shows the maximum 
impact parameter observable with the CFHT MOS FOV, while 
the vertical dotted lines join absorber /galaxy 'pairs' conserva- 
tively identified by their being within 1000 km s _1 of each other 
in redshift. 




Figure 3. Pie diagram showing absorbers and galaxies in the 
LOS of PKS 0122-00. Description of figure as per figure^] 
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Figure 10. Pie diagram showing absorbers and galaxies in the Figure 11. Pie diagram showing absorbers and galaxies in the 

LOS of PKS 1718+481. Description of figure as per figure^] LOS of PKS 2145+06. Description of figure as per figure^] 
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Figure 12. Pie diagram showing absorbers and galaxies in the 
LOS of 3C 454. Description of figure as per figure l2l 



Figure 13. Pic diagram showing absorbers and galaxies in the 
LOS of PG 2302+029 (includes the additional objects from Bowen 
et al. (2001)). Description of figure as per figure |21 
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Figure 14. Pie diagram showing absorbers and galaxies in the Figure 15. Pie diagram showing absorbers and galaxies in the 

LOS of PKS 2340-036. Description of figure as per figure LOS of PKS 2344+09. Description of figure as per figure 
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Figure 16. Pie diagram showing absorbers and galaxies in 
the LOS of LBQS 0107-025A (includes additional objects from 
Rauche, see text). Description of figure as per figure 121 



Figure 17. Pie diagram showing absorbers and galaxies in 
the LOS of LBQS 0107-025B (includes additional objects from 
Rauche, see test). Description of figure as per figure |21 
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Figure 18. Absorption line sample redshift distribution. The 
solid line shows Ly-a lines, while the dotted line shows CIV. 



Figure 19. Galaxy Sample Redshift distribution 



3 ANALYSIS 

3.1 Sample Properties 

The final galaxy sample as described in § EJ] above contains 
685 objects with redshifts. Of these, 49 have 'redshifts' less 
than 500 km s" 1 , and hence are almost certainly stars. As- 
suming this to be true, we can check that our redshift error 
estimates are reasonable. The mean velocity of this sample 
is 5 km s" 1 , with RMS 190 km s~\ The mean of the er- 
ror estimates on the velocities is 230 km s -1 , suggesting we 
might be slightly overestimating our velocity errors, but, to 
be conservative, we make no adjustment for this. 

The absorber sample as described in §EJ] above contains 
815 lines, of which 381 are identified as Ly-a and 54 are 
CIV. Of the CIV sample, there are 25 with unique redshifts 
less than one (i.e. counting doublets as one system). For our 
present analysis we will restrict ourselves to this subsample 
of 406 absorption lines. 

The final redshift distribution of this absorption line 
sample is shown in Figure 1181 The dramatic drop in line 
numbers shortward of z~0.35 is where Ly-a moves below 
the blue edge of the HST FOS G190H grating. In our sub- 
sample, absorption line systems at redshifts less than 0.35 
are generally either metal line systems, or Ly-a detected in 
a GHRS exposure towards LBQS0107-025A or B. 

The final redshift distribution of the galaxies with mea- 
sured redshifts is shown in Figure EH 

In Figure l2771 we show the distribution of the galaxy sam- 
ple with projected distance to the QSO LOS in our adopted 
cosmology. There are relatively few galaxies within 200 kpc 
of the LOS (as would be expected for a uniformly sampled 
distribution). The survey geometry, combined with cosmol- 
ogy, results in a fairly flat distribution out to radial distances 
of ~1 Mpc, with a drop off from there out to 2 Mpc which 
effectively defines the corners of our surveyed volume. 

In FigureEm we show a pie diagram resulting from com- 
bining all of the LOS. This figure should allow the reader to 
understand our survey geometry and the regime over which 
our statistical results apply. In particular one can clearly 
see that our survey covers impact parameters out to 2 Mpc 
for redshifts above 0.45. One can also see that (because the 
KP observations did not include wavelengths shortward of 
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Figure 20. Galaxy Sample distribution in projected distance to 
QSO LOS 



1630A) very few absorber-galaxy pairs are found with red- 
shift less than 0.34. All symbols and lines are as denned for 
Figure EJ 

By comparing the galaxy sample with measured red- 
shifts with the photometric catalogs, we can compute the 
completeness as a function of magnitude for each field in- 
dividually. Because of the fairly small number statistics of 
individual fields, these histograms look fairly noisy, although 
informative. For this paper, we show the merged result from 
combining all of the fields. As can be seen in Figure l2~2l the 
average completeness peaks at around 60% at the bright 
end, falling smoothly to 20% at R~20 and 10% at R~21. 
This completeness can be taken as roughly referring to the 
volume enclosed by the solid line in the bottom panel of 
Figurel2Tl 

3.2 Identifying 'Absorber-Galaxy pairs' 

There are two reasons for trying to identify absorber-galaxy 
'pairs'. One is to investigate the physical relationship (if any) 
between the gas and the galaxies, while the other is to allow 
us to try to tie the velocities of the galaxies to the frame of 
the absorbers. 
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Figure 21. Pie plot showing all Lines of Sight overplotted. De- 
scription of figure as per figure [5] The solid curve in the lower 
panel shows the outer boundary of the survey volume. Also clearly 
visible is the sharp drop in absorber-galaxy pairs which could be 
found produced by the blue edge of the FOS G190H (correspond- 
ing to a redshift of 0.34 for Ly-o ) 



A straightforward way of tying galaxies to absorbers is 
to set limits in impact parameter and velocity difference, 
and to call all galaxies within that distance of an absorber 
'associated'. This begs the question of what to do when there 
is more than one galaxy within this distance of an absorber, 
with common approaches being to just take the closest, as 
defined by some metric combining the velocity difference and 
projected distance, or to accept all such pairs. 

A more complex approach has be en followed by the se- 
ries of papers by Lanzetta and Chen jL anzetta et alj|l99o]. 
19961: IChen et al.l Il998t lOrtiz-Gil et alJ Il999l : IChen et alJ 
2001), who identify 'physical' absorber-galaxy pairs by sep- 
arating them from 'random' and 'correlated' pairs in the fol- 
lowing manner. 'Random' pairs are removed by calculating 
the value of the absorber-galaxy correlation function assum- 
ing a parametric form (kindly provided to us by Hsiao- Wen 



Figure 22. Galaxy sample completeness 

Chen), given the imp act parameter and velocity difference 
jLanzetta et al.l ll99ST) . Pairs are passed to the next stage 
if their calculated correlation amplitude is greater than 1. 
To remove 'correlated' pairs they also then require that the 
impact parameter is less than 200 h" 1 kpc (286 kpc in our 
adopted cosmology) . Given the form of the correlation func- 
tion, and the above additional cut, one can calculate the 
region of impact parameter/velocity difference space within 
which an absorber-galaxy pair would be declared 'physical'. 
We note that this is in practice extremely similar to a simple 
cut at impact parameter of 286 kpc and velocity difference 
of 530 km s" 1 . 

Another route to the above pairing was discussed by 
iMorris et alJ £993), where the three dimensional distance 
between absorbers and galaxies was corrected for the sta- 
tistical fact that objects close in projected distance might 
have relatively large peculiar motions relative to each other 
because they are orbiting in the same potential well. As 
discussed in that paper, it is possible to use the two point 
correlation function to estimate the probability distribution 
of the real three dimensional separation of two objects, and 
then to adopt the expectation value of this probability dis- 
tribution as what was referred to as the 'perturbed Hubble 
flow' distance. 

For simplicity we will only adopt straight cuts in impact 
parameter and velocity difference (shown for example in Fig- 
ure !24H . along with the Lanzetta et al. correlation function 
approach, in this paper. 

3.3 Properties of < 286 kpc. 'Absorber-Galaxy 
pairs' 

In ord er to compare our results with those of IChen et alJ 
( 2001), we have produced pair samples using their approach, 
and also with cuts in impact parameter 286 kpc and veloc- 
ity difference ±530 km s _1 (the region enclosed by dashed 
lines in Figure |2"H . For our sample, these two approaches 
produce identical lists of 13 absorber-galaxy pairs. Given 
our fairly large sample of Ly-cv absorption lines and galaxies 
with redshifts, this small number of 'physical' pairs might 
be seen as surprising, but we re-emphasize that our strategy 
in designing the galaxy redshift survey was not to identify 
the galaxies 'causing' the absorption, but to obtain a sta- 
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Figure 23. Relationship between Equivalent width and impact 
parameter for pairs of absorbers seen in Ly-a and gala xies selected 
as des cribed in the text. Small circles show the data of lChen et alj 
(2001), while larger squares show the new pairs from this paper. 
Also marked are the fiducial EW of 0.3A and impact parameter 
of 257 kpc (see text) 



tistically well-understood sample of galaxies over a region 
large in impact parameter. For this reason we did not try to 
concentrate our efforts at small impact parameters. Indeed, 
this would have been a rather inefficient use of telescope 
tim e and MPS spectro graph real-estate. 

IChen et alJ (1200 ill use their sample of 34 galaxy- 
absorber pairs to conclude that L* galaxies have tenuous 
gaseous halos of column density ^10 14 cm" 2 , with radius 
180 h _1 kpc and covering factor near unity, and that the 
halo radius scales with galaxy B-band luminosity as Lg 39 . 
Their Figures 3 and 4 make the case for this conclusion. In 
Figure 1231 we show a reproduction of their Figure 3 (con- 
verted to our adopted cosmology) along with the additional 
pairs from our sample. Also marked are the fiducial EW of 
0.3A (corresponding roughly to a column density of 10 14 
cm~ 2 ) and the L* impact parameter of 257 kpc (180 h _1 
kpc with h=0.7). 

In response to an interesting question from the referee 
of this paper, we have investigated the sub-sample of pairs 
found with large rest EW (taken to mean larger than lA) 
and large impact parameter (taken to mean larger than 100 
kpc). There are 5 such pairs in our sample. With such a 
small sample, it is impossible to draw any statistically con- 
vincing conclusions, but it is notable that all 5 galaxies so 
identified have at least one other galaxy close to it in red- 
shift. The clearest example is the absorber galaxy pair in the 
PG 2302+029 LOS at redshift 0.59 where there is a clear 
group or wall of galaxies at this redshift (see Figure "TBI . 
One possible inference, which could only be confirmed with 
a much larger sample, would be that large EW, large im- 
pact parameter 'pairs' are identified in somewhat overdense 
regions. 

The additional pairs from our sample do significantly 
increase the number of large column density pairs with sep- 
arations larger than 100 kpc, but the incomplete nature 
of our sample (and indeed of most samples of this nature) 



does leave open the possibility that a galaxy of comparable 
brightness might be found closer to the LOS. 

We also looked at the rest frame B-band luminosities 
of the galaxies that are paired in this manner, finding that 
these were not significantly different from the luminosities of 
all the galaxies found in pairs with separation up to ±5000 
km s _1 and impact parameter up to 2 Mpc. I.e. there is 
no evidence from our sample that more luminous galaxies 
are more likely to be found in absorber-galaxy pairs. How- 
ever, the large uncertainties in the derived absolute B-band 
magnitudes we tabulate, that results largely from an uncer- 
tain k-correction, mean that any such trend, even if present, 
could be undetectable using the data in this paper. 

As a result of these inconclusive results, we have cho- 
sen to abandon the approach of trying to identify a single 
'physical' galaxy-absorber pair, and in the next section will 
concentrate on the statistical properties of all the galaxy- 
absorber pairs within a fairly generous window. 

3.4 Properties of All Absorber Galaxy Pairs 

with velocity difference less than 1000 km s _1 

We can now analyse the statistical properties of the relative 
distributions of the gas and galaxies along the 16 lines of 
sight. The basic sample we will use will be taken from all 
possible pairs of absorbers and galaxies with velocity differ- 
ences less than +/- 5000 km s" 1 . For the CIV systems we 
have only taken the stronger line of the doublet. We note 
that therefore individual absorbers or galaxies can appear 
multiple time as different pairs. The entire sample is shown 
in Figure 124*1 This complicated plot shows both the regions 
we have determined to contain an excess of absorber-galaxy 
'pairs', and also a region beyond +/- 1000 km s _1 which 
we will use to determine the 'background' value for random 
pairings. Our sample is large enough to break into bins of 
impact parameter (shown divided by horizontal lines in the 
figure) . We also measure absorber-galaxy pair over-densities 
for three different velocity ranges +/- 200, 500 and 1000 km 
s _1 . We show the illustrative box at impact parameter of 
286 kpc and velocity differenc e of 530 km s~ which a pprox - 
imates the r egion s elected by [Lanzetta^t aD (^9^_[w9^: 
IChen et alJ (Il998l) : lOrtiz-Gil et all <ll999l) and IChen et*al] 
fcOOlli as containing physical pairs. 

In Figure 1251 we show the above data projected onto 
the x-axis and displayed as a histogram. We also show the 
mean background values for the Ly-a and CIV samples. Two 
thing can be readily seen from this figure. First that there is 
a strong over-density of pairswith a small velocity difference 
in both Ly-a and CIV, and secondly that the velocity differ- 
ences seen in CIV are somewhat smaller than those seen in 
the Ly-a . The RMS of the velocity differences between all 
absorber-galaxy pairs seen in Ly-a with velocity difference 
between ±1000 km s -1 is 440 km s _1 (191 pairs), while the 
RMS for absorber-galaxy pairs seen in CIV with velocity 
difference between ±1000 km s _1 is 350 km s _1 (31 pairs). 

We also show in Figure 1261 how the pairs in the two 
other plots are distributed in redshift. As a result of the 
wavelength coverage of the FOS spectra providing our ab- 
sorption line sample, there is a clear break in the sample at 
z=0.29. For the following analysis we will consider both the 
whole sample, and also just the pairs with z>0.29 to check 
whether the results are being affected by the low redshift 
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Figure 24. Distribution of velocity differences and impact pa- 
rameter between absorbers seen in Ly-a or CIV and galaxies se- 
lected as described in the text. Absorbers seen in Ly-c? only are 
marked by open circles, absorbers seen in CIV are marked as filled 
circles. 



Figure 26. Distribution of redshift and impact parameter for 
absorber-galaxy pairs. Pairs with absorbers seen in Ly-a only are 
marked by open circles, pairs with absorbers seen in CIV are 
marked as filled triangles. 
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Figure 25. Distribution of velocity differences between absorbers 
seen in Ly-a or CIV with galaxies selected as described in the text. 
The upper line shows the pairs with absorbers seen in Ly-a only 
while the lower line shows pairs with absorbers seen in CIV. 
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Figure 27. Measured Excess of Absorber-Galaxy Pairs relative 
to the Background 



Absorber-Galaxy Excess z>0.29 



sub-sample (which mostly comes from the LBQS 0107-025A 
and B LOS). 

Given this sample, we can now pose the question: How 
large is the absorber-galaxy pair over-density compared to 
the background, and how significant is the detection of this 
over-density? In order to measure the latter, we consider the 
numbers of pairs seen in the various bins in impact param- 
eter and velocity difference in Figure 121 to be drawn from 
a Poisson distribution with an expected mean given by the 
'background' density measured at the same range of impact 
parameter but with velocity differences between -5000 and 
-1000 km s" 1 and between +1000 and +5000 km s" 1 . The 
results of the analysis are shown in Figures 1271 1281 1291 and 

no 

In Figures 1271 and 1281 we show that (as expected by a 
visual inspection of Figure 1241 there is a substantial over- 
density of absorber-galaxy pairs at impact parameters less 
than 500 kpc. The fractional over-density is larger for the 
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Figure 28. Measured Excess of Absorber-Galaxy Pairs with 
z>0.29 relative to the Background 
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Figure 29. Significance level of the Measured Excess of 
Absorber-Galaxy Pairs relative to the Background. The number 
plotted on the y-axis is the probability that the number of pairs 
found within velocity range given in the sidebar would be found 
given the mean background and no real excess. 
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Figure 30. Significance level of the Measured Excess of 
Absorber-Galaxy Pairs with z>0.29 relative to the Background. 
The number plotted on the y-axis is the probability that the num- 
ber of pairs found within velocity range given in the sidebar would 
be found given the mean background and no real excess. 



systems seen in CIV than for the Ly-a systems. At larger 
impact parameters, the fractional over-density is less. It is 
reasonable that the strength of the over-density generally 
drops as one increases the velocity difference range within 
which one counts pairs. 

In Figures l2§l and l3Tjl we show that the over-densities no- 
table in Figures ITFl and l^Hl have a <1% probability of arising 
by chance out to impact parameters of at least 1. 5 Mpc. This 
is con sistent, for example, with the models of iDave et alJ 
(1999). Figure 18 in that paper for example shows corre- 
lations in absorber-galaxy properties out to these sorts of 
scales. 

Note that the probabilities reported assume that the 
expected (or background) value is known perfectly. In prac- 
tice there is a (Poisson) error on our measurement of this 
background, which could be folded in to the analysis, but 
has not been in order to avoid overly complicated plots. 



Figure 31. Distribution of velocity difference and impact pa- 
rameter between galaxies and galaxies selected as described in 
the text. 



As noted in § 12.21 the identifications of the lines in the 
PG 1718+481 LOS are not complete. We have rerun the 
above statistical tests including all the unidentified lines as 
Ly-a and find that this makes no significant difference to 
the results reported in this section. 

We also remind the reader that the FOS spectral reso- 
lution of 230 km s _1 means that no pairs of absorption fea- 
tures can be identified with velocity difference smaller than 
this. The low density of Ly-a absorbers at these redshifts 
make the likelihood of two such lines coinciding by chance 
low, although clustering can obviously chan ge this conclu- 
sion. T his issue is discussed in more detail in ljannuzi et al.l 
(1998). We also note that the galaxy redshift uncertainties 
(also of 230 km s _1 ) will artificially move galaxies across 
the boundaries of any velocity difference cut. The general 
trend will be (a) to weaken any measured correlation, and 
(b) for this weakening to be stronger for the smaller velocity 
difference bins. 

An obvious question is whether the observed results are 
simply due to the absorbing gas being in the halo of a galaxy, 
and the detected pair excess being due to the galaxy-galaxy 
correlation function. This is tested using the galaxy sample 
and forming a similar sample of galaxy-galaxy pairs to the 
absorber-galaxy pairs. In order to approximately match the 
absorber galaxy sample, we choose galaxies from each field 
within 500 kpc of the QSO LOS, and collect up all resulting 
pairs with that sample. This subsample should then be able 
to find pairs in the same volume of space as the absorber 
sample. Figures UTTl and l3l?l are the equivalent of Figures 
and 1261 but for the galaxy-galaxy pairs. The over-densities 
found for galaxy-galaxy pairs were plotted in Figures and 

The above analysis shows that the galaxy-galaxy pair 
over-density is significantly higher than the absorber-galaxy 
pair over-density seen for the Ly-a systems, but compara- 
ble to that seen for the absorber systems seen in CIV. This 
result could have several interpretations, but one might be 
that the absorbing gas seen in CIV (and which has obviously 
been polluted by outflows from galaxies) is still in close prox- 
imity with galaxies, while systems seen in Ly-a are more 
widely distributed, albeit still correlated with the galaxy 
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Figure 32. Distribution of redshift and impact parameter be- 
tween galaxies and galaxies selected as described in the text. 

distribution. 5 Another possible interpretation is that we are 
just seeing the column density correlation with impact pa- 
ra meter noted by ot her authors and modelled for example 
bv lDave et alJ (ll99STl . Separating hydrogen column density 
from metallicity as the underlying driver of these correla- 
tions needs further work and will benefit from higher resolu- 
tion UV spectroscopy (i.e. data from the HST archive taken 
with the Space Telescope Imaging Spectrograph (STIS)). 

3.5 The link to simulations and theory 

One of the many benefits provided by the powerful cosmo- 
logical simulations begun in the 1990s and continuing into 
the present (see introduction for list of references), is the 
ability to simulate the evolution of the physical properties 
of Ly-a absorbers as a function of redshift and cosmic time. 
This helps tremendously when trying to make use of diverse 
observational data sets in order to understand what physi- 
cal processes influence the physical properties and evolution 
of the IGM and galaxies. The theoretical models aid us in 
identifying the high redshift progenitors of structures ob- 
serv ed at low redshift. For example, i t is from such models 
fe.g. lDave et alJ lll999l) : ISchavel d200ll) ') that we understand 
that a given relative over-density in the (dark) matter dis- 
tribution will be associated with a different column density 

5 When this paper was near to submission, an interesting short 
letter appeared on astro-ph bv lChen et"ai] <2005l) describing a sin- 
gle LOS with high quality STIS UV spectroscopy. They perform 
a correlation analysis with a sample of 61 galaxies within 1 h _1 
Mpc at redshifts less than 0.5. Some of their conclusions will need 
confirmation with a greater number of LOS and a larger galaxy 
sample, but their conclusion that the absorber-galaxy correlation 
is entirely due to the emission line (i.e. star forming) galaxies is 
potentially very exciting. If confirmed, then this would support a 
seductively simplistic model with star forming galaxies and high 
column density absorbing gas co-existing in the filamentary struc- 
ture now so familiar from the simulations, while the more highly 
clustered absorption line galaxies lie in the knots and intersec- 
tions where the 10 4 K gas needed for detection in Ly-a is absent 
or destroyed (possibly by conversion into the much sought after 
WHIM). Our data set can be used to look for confirmation of this 
picture. 



of neutral hydrogen as a function of redshift. Ly-a absorbers 
with Nhi approximately 10 14 cm -2 at z = are in regions of 
over-density that at z = 3 would have neutral gas columns 
closer to 10 16 cm~ 2 JPave et alj H999L We note that this 
does not imply that regions of neutral column density 10 16 
cm -2 at z—S evolve into regions with neutral column den- 
sity 10 14 cm -2 at z—0, but merely that these two different 
column densities of neutral material flag regions at those 
two epochs with the same dark matter overdensity. 

A natural next step for this work will be to compare 
the statistics we have derived from our observations with 
statistics measured as close as possibly identically on sim- 
ulations which include feedback, with a goal of verifying or 
ruling out the feedback assumptions they have made. This 
is a challenging project, as the simulation must faithfully in- 
clude both low density inter-galactic gas, and also accurately 
predict the locations and physical properties (stellar lumi- 
nosities and star formation histories) of the highly-collapsed 
galaxies. The redshift range from z=l to the present day 
is well known to be hard to model at high resolution, be- 
cause of the large dynamic range needed. That said, several 
of the references to modelling work given in the introduction 
suggest that this problem is tractable. 



4 CONCLUSIONS 

A goal of our work is an improved understanding of how gas 
is transformed into galaxies, and how those galaxies in turn 
influence the gas around them (e.g. by winds, jets or radi- 
ation). Studies of absorbers and galaxies at high redshift 
have already yielded evidence for such mechanisms play- 
ing a role in the f ormat io n and evolu t ion of galaxies (e.g. 
lAdelbereer et alJ lEoo3): ICroft et al.1 (120021); iFrank et alJ 
j2003h : lBruscoli et alJ J2003ft and iKollmeier et al.l (l2003l) L 

We have focussed on the second half of the history of the 
Universe, i.e. at redshifts less than one. This is a period in 
which the star formation density of the Universe is thought 
to be dropping away from its peak, the rapid evolution of the 
properties of the Ly-a absorbers has slowed, and when the 
Hubble sequence of galaxies is well established. Despite this 
apparent middle-aged placidity, there should remain clear 
evidence of the wild excesses of youth, and that is what 
we claim to have measured. An advantage of focussing on 
this (large) fraction of the history of the Universe is that 
the galaxy population can be studied in detail, and with 
confidence that a representative sample has been obtained. 

There is evidence from the KP data set alone for cluster- 
ing of some fraction of low redshift (z < 1) Ly-a absorbers 
on velocity scales < 300 km/sec. This evidence is that ap- 
proximately 10% of the KP Ly-q absor bers are "resolved" 
in the FOS spectra jjannuzi et alj|l993) . 

At higher column densities, the Lyman limit systems 
contained in the KP sample all have extensive associated 
metal line systems. For all the extensive metal line systems 
(multiple absorption lines from strong resonance line pro- 
ducing species, e.g. Si, O, N, etc.) for which we could check 
for the presence of an associated Lyman Limit system, we 
detect such a system. In other words, the absorption line 
systems that are unambiguously higher column density sys- 
tems (at low redshift) are unquestionably more like the gas 
producing the ISM of a galaxy than like the gas one might 
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expect in a void. These same absorbers also often have asso- 
ciated OVI, or NV, or other evidence of more highly ionized 
components i.e. multi-phase gas again, like a galaxy. 

Our conclusions for the more moderate column density 
systems that dominate our current sample can be summa- 
rized as: 

• A correlation between absorbers and galaxies has been 
detected out to impact parameters of at least 1.5 Mpc. 

• The strength of the absorber-galaxy correlation is 
weaker than the correlation between galaxies and galaxies. 

• The velocity differences seen between galaxies and ab- 
sorbers with detected CIV are typically smaller than the ve- 
locity differences seen between galaxies and absorbers seen 
in Ly-a only. 

• The above is consistent with absorbers being a mixed 
population containing some pristine 6 material probably in- 
falling for the first time along the filaments predicted by 
current models, and some contaminated material produced 
by outflows which in general lies closer to galaxies in veloc- 
ity. 

• The above qualitative picture needs to be fleshed out 
by comparing the numerical strengths of the observed cor- 
relations and relationships with SPH or AMR modelling. 

While we have so far only made a fairly superficial anal- 
ysis of the (no-doubt complex) relationship between gas and 
galaxies at redshifts less than one, we feel the above results 
are very encouraging, and that they suggest that there is 
considerable room for further observation and modelling of 
this crucial interaction. 
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